Non-Fermi liquid and topological states with strong spin-orbit coupling 
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We argue that a class of strongly spin-orbit coupled materials, including some pyrochlore iridates 
and the inverted band gap semiconductor HgTe, may be described by a minimal model consist- 
ing of the Luttinger Hamiltonian supplemented by Coulomb interactions, a problem studied by 
Abrikosov and collaborators. It contains two-fold degenerate conduction and valence bands touch- 
ing quadratically at the zone center. Using modern renormalization group methods, we update and 
extend Abrikosov's classic work and show that interactions induce a quantum critical non-Fermi 
liquid phase, stable provided time-reversal and cubic symmetries are maintained. We determine 
the universal power-law exponents describing various observables in this "Luttinger Abrikosov Be- 
neslavskii" state, which include conductivity, specific heat, non-linear susceptibility and magnetic 
Gruneisen number. Furthermore, we determine the phase diagram in the presence of cubic and/or 
time-reversal symmetry breaking perturbations, which includes topological insulator and Weyl semi- 
metal phases. Many of these phases possess an extraordinarily large anomalous Hall effect, with the 
Hall conductivity scaling sub-JinearJy with magnetization, a^y ~ M"'^^. 



Divining the nature of critical non-Fermi liquid (NFL) 
phases of electrons is one of the most outstanding prob- 
lems in correlated electron systems, such as cuprates, 
pnictides, and heavy fermion materials [IHl]- Theory, 
spurred largely by NFL behavior in the cuprates, has fo- 
cused on models with large Fermi surfaces, appropriate to 
these materials. Though there have been promising tech- 
nical advances [S], this approach has not provided clear 
resolution of experimental puzzles. Recent theory and 
experiment have uncovered a new frontier for correlated 
phenomena: Mott correlation physics in materials with 
strong spin-orbit coupling (SOC) [Si [7]. Of particular 
interest in this regard are the 5d transition metal ox- 
ides, where compelling evidence has been built for Mott 
phenomena driven cooperatively by SOC and Coulomb 
interactions [3HS1- Here we consider the possibility of 
distinct NFL states in the strong SOC regime. 

We are particularly motivated by the pyrochlore iri- 
dates, A2lr207, where A is a lanthanide element llOj. 
These materials display T > metal-insulator transi- 
tions, with the metal-insulator transition temperature 
decreasing with increasing A radius. The progression 
culminates with Pr2lr2 07, which is a highly unconven- 
tional metal down to the lowest temperatures [TTl [T^] . It 
displays logarithmic NFL behavior of the magnetic sus- 
ceptibility, and a remarkable enormous zero field anoma- 
lous Hall effect ( AHE) , in the absence of any measurable 
magnetization [TTl [H] ■ 

With these studies as motivation, we utilize prior stud- 
ies of the electronic structure of pyrochlore iridates [TBl- 
[T5] to show that a minimal description for the electronic 
states in their paramagnetic phase is a storied Hamilto- 
nian from semiconductor physics: the Luttinger model 
of inverted band gap semiconductors [16]. This model 
has gained recent notoriety for its relevance to HgTe, the 



starting material for some topological insulators [T7U20j . 
While HgTe is a weakly correlated material where band 
structure alone provides a sufficient description, in the 5d 
materials the Luttinger Hamiltonian must be supplanted 
by interactions. 

In this paper, we carry out such an analysis, rediscov- 
ering and extending a storied analysis by Abrikosov and 
Beneslavskii of Coulomb forces on the Luttinger prob- 
lem [in [22] ■ Using modern renormalization group (RG) 
techniques, we confirm Abrikosov's conclusion that long- 
range Coulomb interactions convert the quadratic band 
touching into a quantum critical NFL, prove the stabil- 
ity of the state within an e expansion, and calculate the 
full set of anomalous dimensions characterizing the state. 
Consequently, we call the resulting phase a Luttinger- 
Abrikosov-Beneslavskii (LAB) state. While the LAB 
phase is stable in the presence of time-reversal and cubic 
symmetries, we show that it is a "parent" state for other 
exotic states that can be reached by breaking either or 
both of these: metallic and double- Weyl semi-metallic 
phases with enormous anomalous Hall effects (AHEs), 
and topological insulators. We discuss the implications 
for the iridates at the end of the paper. 

The LAB phase itself has striking properties. Its NFL 
nature is revealed directly by algebraic singularities in 
the electron spectral function (probed in angle-resolved 
photoemission) and in optical conductivity, as well as 
indirectly through many power law thermodynamic and 
response functions. 

We now turn to the exposition of these results. We 
consider the paramagnetic band structure based on prior 
work [T3HIS], which argued that the states at the zone 
center (F point) near the Fermi energy are comprised 
of the four-dimensional representation, which can be de- 
scribed by "angular momentum" operators J (which are 
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FIG. 1. Phase diagram of the perturbed LAB in the space 
of renormaUzed strain to Zeeman field ratio, A = {S/H)u, 
versus cubic Zeeman angle, 0. Here"Weyl" denotes the (dou- 
ble) Weyl semimetal, "Ins." insulator, and "Metal" a metallic 
phase which has Weyl points shifted from the Fermi energy 
in the region below the dashed line. For H — 0, the insulator 
is a topological insulator. 



j = 3/2 matrices) transforming as the T2 representation 
of the cubic group. In our minimal model, we assume 
only these states close to F are important. Then k ■ p 
theory and cubic symmetry determines the band struc- 
ture in its vicinity to be precisely described by the Lut- 
tinger Hamiltonian with three effective mass parameters 
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This describes doubly degenerate bands with energy 
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where Pc{k) — 
(4M,Mo)/(4M,-5Mo; 



K - E.^, fc?fc| and Mo = 
Henceforth we assume Mq > m, 
which describes conduction and valence bands touching 
quadratically at E = k = 0, where the chemical potential 
for the undoped material crosses. 

The LAB is obtained by adding to this the long-range 
Coulomb interaction. We implement the latter by a 
scalar potential ip, which in the Euclidean path integral 
formalism gives the action 
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with Ho = ^o(— «V) and cq — l/An. Here i/' is a four- 
component spinor, but subsequently we will artificially 
add an additional U{Nf) flavor index, which allows a 
check on our calculations by large Nf methods; the phys- 
ical case is Nf = 1. Eq. ^ contains in addition to the 



three mass parameters, the Coulomb coupling constant 
e. For e = 0, scale invariance is manifest, with the scaling 
dimensions [x~^] = 1 , [r^i] = z , [1/^] = f , [^] = z - 2, 
[if] = {d + z — 2)/2. Here we introduce the dynamic crit- 
ical exponent (z), which is naturally z = 2 with e = 0, 
but will become non-trivial with interactions. 

Directly in the physical case d = 3, the dimension of 
the coupling constant is [e^] — 1, so Coulomb interactions 
are strongly relevant. Therefore we employ the e = A — d 
expansion to control the RG analysis. As familiar from 
quantum electrodynamics, three one loop Feynmann di- 
agrams contribute to leading order in e: the fermion self- 
energy, boson self-energy, and vertex correction. Here we 
show that the relevance of Coulomb interactions signals, 
rather than a flow to strong coupling and a symmetry 
breaking instability, the formation of a new stable inter- 
acting flxcd point, which describes the critical non- Fermi 
liquid LAB state (Abrikosov's analysis tacitly assumes 
this stability). 

The RG is carried out perturbatively in e, but non- 
perturbatively in the mass parameters. Thus a full treat- 
ment gives non-trivial and complete beta functions for 
the two dimensionless mass ratios m/Mo, m/Mc] these 
are given in the Supplementary Material. The analysis of 
the full RG shows, however, that there is a single stable 
isotropic fixed point corresponding to tji/Mq — m/Mc = 
0, so for simplicity we quote in the main text only the 
results in the vicinity of this point. 

In this limit, the leading contribution to the bosonic 
self energy becomes 
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where we took the w — >■ limit because frequency de- 
pendence is subdominant. The divergence should be ab- 
Pc{k), (1) sorbed by rescaling the bosonic field, (p — t- e^'"'''^(^ upon 
reduction of the hard momentum cutoff A — >■ e~'*^A, 
which defines the RG parameter £. This gives the bosonic 
anomalous dimension rji, = 2NfU [241, where the dimen- 
sionless coupling constant is u 



a 2^\i-d , which has 
the physical meaning in d = 3 of the ratio of the real 
space cutoff to the effective Bohr radius. The frequency 
dependence of the one loop fcrmionic self-energy and the 
vertex correction both vanish, the result of a Ward iden- 
tity. For fc 7^ 0, the fermion self-energy gives mass cor- 
rections, e.g. 5{l/m) = 8u/(15m) x d£ to leading order. 
Detailed analysis is given in the Supplementary Material. 

Given these calculations, we choose z = 2 — 8u/15 to 
keep the mass m fixed, which gives the RG equations, to 
lowest order in m/Mcm/Mf): 
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From the first equation above we find tlie fixed point 
coupling and hence dynamical exponent, 

u* — ^ e , z = 2 e. (5) 

307V/ + 8 ' 15iV/+4 ' ^ ' 

and since u* > 0, the second line in Eq. Q implies both 
itl/Mq and m/Mc are irrelevant. This establishes the 
existence and nature of the stable, isotropic fixed point 
describing the LAB phase. As a check, we have carried 
out a large Nj expansion, which gives the same bosonic 
anomalous dimension as in the e expansion at the one- 
loop level, supporting the stability of the LAB phase. 

The presence of the stable interacting fixed point can 
be understood physically as a balance of partial dynam- 
ical screening of the Coulomb interactions by electron- 
hole pairs and mass enhancement of the same quasipar- 
ticles by pairs. This situation is in sharp contrast to the 
case of a vanishing indirect band gap, for which to lead- 
ing order in the long-range Coulomb interaction electrons 
and holes are separately conserved, so there is no screen- 
ing by virtual electron-hole pairs, and exciton formation 
destabilizes the putative gapless state[25]. 

Using the RG, we can evaluate the anomalous di- 
mension of any physical operator. By charge conser- 
vation, [V'^V'] — d. Due to the isotropy of the fixed 
point, there are only two non-trivial values for the 
other charge-conserving fermion bilinears. We obtain , 
[ip'^Taip] =d + r]i , [ip^Tab'ip] = d + r]i2, , where are the 
(time-reversal invariant) Dirac gamma matrices. Tab = 
— |[Fo,F{,] are time-reversal odd, and a, & = 1, 2, • • • , 5. 
Using the standard operator insertion technique, we find 

^1 ^ ~ i5Nf+i ^ ^^'^^ '^12 = ~ i5Nf+4: ^- These operators 
describe many physical observables, e.g. the "angular 
momentum" operator ~ {—T34—^Ti2)^- The nega- 
tive anomalous dimension of these operators suggests the 
a schematic picture of power-law excitons due to electron- 
hole attraction. For pairing channels, we find positive 
anomalous dimensions, consistent with this view. The 
local pairing channel has -qpairing = ^ = 30Af%8 ^- 

Using these results, we obtain thermodynamic re- 
sponses such as the specific heat, Cy ~ T'^/'^ sa T^-^ 
and the spin susceptibility x{T) ~ a -f 5T(rf-^+2'7i2)/2: ^ 
a + hT^'^, with some constants a, fe. Interestingly, 
the non- linear susceptibility X3 — M / dH^\j^_^ ^ 
rp-(3z-4r}i2-d)/z ^ -J,-!.? (jiyerges, as in spin glasses but 
with completely different physics. Comparing the scaling 
of current and electric field gives the usual result [aij] = 
d — 2. Consequently, the temperature and frequence de- 
pendence of the conductivity is a 

and a clean, undoped LAB is therefore a power-law insu- 
lator. 

We now turn to the effect of applied strain and Zee- 
man field upon the LAB. These perturbations break 
cubic/time-reversal symmetries, and thus destabilize the 
LAB. Due to the isotropic nature of the LAB fixed point. 



the response to the Zeeman field alone is to leading or- 
der independent of its direction (the cubic mass 1/Mc 
can be "dangerously irrelevant", however - see below), 
so we take it to lie along the (001) direction. We con- 
sider for simplicity tetragonal strain which preserves C4 
rotation about this axis (in the absence of Zeeman field, 
the direction of strain is again unimportant). This leads 
to the perturbations 

%' = -5{Jl -\)- H{cos{e)J, + sm{9)J^), (6) 

where S parametrizes the strain, h is the Zeeman field, 
and 9 controls the strength of the cubic Zeeman term 
allowed by the cubic symmetry [551 HZ] • Using the RG 
results, the dimensions of these perturbations are [6] = 
z — rji « 2.1 and [H] = z — 7712 « 1.9; i.e. strain is slightly 
enhanced while Zeeman field is slightly suppressed by 
interactions. However, both dimensions are positive and 
close to 2, so that they are strongly relevant. They fiow 
to strong coupling under the RG, and the fate of the 
system must be re-analyzed in the limit. 

To do so, we assume, and check self-consistently, that 
interactions have weak effects at strong coupling, and 
simply solve the quadratic Hamiltonian (with tu/Mq = 
m/Mc = 0) in the presence of the renormalized T-l' . 
The result depends upon the dimensionless quantities 
9 and the renormalized coupling ratio A = {S/H)ji ^ 
j/7j(^-m)/(2=-r,i2). For 7J = (A = 00), we have time- 
reversal invariance, and we recover the known result that 
strain S > induces a gapped, 3d TI phase, as observed 
in HgTe [19]. The situation in applied Zeeman field is 
more interesting. Notice that for k = kz, Jz is a good 
quantum number, and there is no level repulsion between 
bands of different . This allows (non-degenerate) bands 
to cross along this axis, which indeed occurs when |A| is 
not too large. Further analysis in the Supplementary Ma- 
terial shows that these crossings correspond to a pair of 
double Weyl points, with linear dispersion along the z axis 
and quadratic dispersion normal to it. These points are 
strength ±2 monopoles in momentum space. Away from 
the kz axis, electron and hole pockets may accidentally 
cross the Fermi energy. If this does not occur, one has 
a pristine double Weyl semimetal, which occurs for the 
angular range 9i < 9 < 92, where 9i — — tan"-'- ( rf^_^ ) 

and 6*2 = tan~-'-( ^^^'^ ) for (5 = 0, as shown in horizontal 
axis ofFig.[lj WhenO < |A| < 00, we observe insulating, 
double Weyl semimetal, and Weyl metal (with coexisting 
electron-hole pockets) phases, as shown in Figure[l] Note 
that in all these phases, the Coulomb interactions become 
either unimportant (in the insulator), screened (in the 
metal) , or marginally irrelevant (in the Weyl semimetal) , 
justifying our treatment of the phase diagram to a first 
approximation. 

More subtle effects may make small modifications to 
this picture. Coulomb interactions can destabilize some 
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of the quantum phase transitions in Figure [T] leading to 
intermediate phases. When the magnetic field is applied 
along a low symmetry axis, the double Weyl points can 
split into multiple single Weyl points, once the effects of 
the cubic mass 1/Mc is included, which is dangerously 
irrelevant in this case. 

A striking experimental consequence of this phase di- 
agram is the AHE due to the Zeeman field, which could 
originate either from an external magnetic field or as an 
exchange field due to local moments in the material. The 
latter is particularly interesting in light of the experi- 
mental results on Pr2lr2 07, which shows a large AHE 
in a regime where the magnetization M is immeasur- 
ably small [m [H]. On symmetry grounds, a^y ^ 
implies ^ 0, but evidently tr^j, is unusually large 
relative to M . This behavior is in fact characteristic of 
the LAB: since the Hall conductivity has dimensions of 
inverse length, we expect Uxy ~ H'^/i'^-v^-^) « H^-^^, In 
the situation relevant for Pr2lr2 07, the Zeeman field is 
generated by Kondo exchange with the Pr moments, so 
H ~ JrM, with M the (dimensionless) Pr magnetiza- 
tion, which implies a highly unconventional sublinear de- 
pendence of axy on M for the pristine LAB. If the Fermi 
level is displaced from the band touching by an amount 
ep, then we expect, treating the above power as a square 

root, a^y ~ ^^J^^j^S^H/ep), where 5 is a scahng func- 
tion (see Supp. Mat.). This gives an order of magnitude 
quantitative estimate 




JkM > ep 
JkM < ep 



(7) 

where Ry=13.6eV is the Rydberg, and rrie is the electron 
mass. The lower regime gives axy ^ O.lfl^^cm^^ , within 
an order of magnitude of observations in Pr2lr207,^2j for 
parameters m ~ 20me (estimated from the calculations 
in Ref. lTS)) , ep ~10meV, Jk ~100K (estimated from the 
measured Curie- Weiss temperature), and M ~ 0.01. 

Another interesting experimental observation in 
Pr2lr207 is a diverging magnetic Gruneisen number 
Th = ^ §jj I g in the zero field limit pSMTj . In a purely 
electronic system with no local moment contribution to 
the entropy, we can readily obtain the behavior of the 
LAB in the low temperature limit. 



r.(i/,T).-4^1, 

yo(z + rjHj H 



(8) 



which depends upon the exponent ?/o defining the tem- 
perature dependence of the specific heat, C ^ T^°, of the 
LAB phase. In the Weyl metal, j/o = Ij and Th < 0, 
while for isolated double Weyl points, yo — 2, and 
> 0. Thus one may imagine a sign change of the 
Gruneisen number when field or strain is varied. Note 
that in Pr2lr207, there is certainly a large local moment 



contribution to the entropy, so that Eq. ([8]) is not liter- 
ally applicable. Nevertheless, the LAB physics may play 
some role in this quantity. 

In conclusion, we have described a novel NFL phase 
occurring in correlated strong SOC systems, with a natu- 
ral connection to the pyrochlore iridates. Even with weak 
correlation effects, some of the phenomena discussed here 
can be observed with only minor modifications, and it 
would be interesting to search for them in HgTe. Fu- 
ture theoretical studies should include more comprehen- 
sive treatment of breaking of cubic symmetries, and the 
effects of disorder and doping. 
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SUPPLEMENTARY MATERIAL 
Representations of the Hamiltonian 

In this section, we provide more information about 
the different representations used in the main text. The 
Hamiltonian reads 

Hoik) = aifc2 + a2{k ■ jf + ^3(^2 + fc^ j2 + fc^ j|) 



e , (fP-fc.J)2 ^klJl + kpl + klJl) 
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The first line uses the conventional Luttinger parameters 
(0^1,2.3) in the j = 3/2 matrix representation, and the 
second line is the form used in the main text. For the 
purpose of computations, it is convenient to introduce 
the Clifford gamma matrices (F^) in the third line as in 
the paper by Murakami et al. |23) . 
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di{k) = -Vikykz , d2{k) = -VSk^k^ , dsik) = -Vsk^ky, 
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d,{k) = —{2kl-kl-kl). (10) 



It is straightforward to relate the masses used in the main 
text and the Luttinger ai parameters. This can be done 
by expressing the spin operators in terms of gamma ma- 
trices, using the equalities 
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Weyl semimetal 



Here we consider how the Weyl semimetal appears, 
taking for simplicity the case 5 = 0,6i < 9 < 92- It is 
straightforward to generalize this to the cases with 5 7^ 0. 
For convenience, we set m = 1/2. For k = (0,0, kz), the 
Hamiltonian in the presence of the Zeeman field, h (along 
the same axis), becomes 

nikzz) = kl{^ - Jl) - H{cos{9)Jz + sin(6') Jf). 

Clearly, the energy eigenstates can be labeled by both 
kz and Jz — ±1/2, ±3/2. One can readily see that level 
crossings occur between the two bands with = 1/2 and 
Jz = —3/2. In the vicinity of these crossing points, the 
other states with Jz = —1/2, +3/2 can be discarded, and 
the Hamiltonian projected onto the two level subspace 
of low energy states. We introduce Pauli matrices in 
this subspace, so that = [^(^l — | — |)(— 1|, = 
(r- + ^Ty)/2 = |i)(-f| = (r-)t, and = + | - 

|) ( — 1 1 , which is the identity matrix in the 2x2 subspace. 
We define two energy parameters, e_3/2 = -ff(| cos{9) + 
f sin(6')) and ei/2 = -H{\ cos{9) + | sin(6l)). Then the 
reduced Hamiltonian becomes 

= e+r" + (e_ + dslfc))^" + d4(A)T" - d3(^)^^ 
e± = (e-3/2±ei/2)/2. (12) 

There are level crossing points at k^ — ky ^ and kz — 
±K, with K = y/eZ. We expand around these points, 
letting kx = Px, ky = py and kz = ±K + pz , which gives 
'H2{±Kz +p) = e+T^ + T~L2{p), with, to leading order in 
P, 

= TvpzT' + d^{p)T- + d^{p)Ty , 
with V = 2yfeZ. The energy spectrum is 



e{p) = ±\Iv-^pI + -{pI 
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(11) 



We see that the electrons disperse linearly along the 
field and quadratically orthogonal to it, near the touch- 
ing point. This can be understood as a consequence 
of 4-fold rotational symmetry around the z axis. Since 
\Jz) — ?> e"^'^^/2| Jz) under such a rotation, the operators 
carry a net angular momentum of ±2, and therefore 
must couple to the "d-wave" combinations of p^. andpj,, 
which are precisely given by d^(p) and d^ij)). 

Though the quadratic dispersion normal to the field 
is due to symmetry, the touching itself has a topological 
character. To see it, it is convenient to define the reduced 
Hamiltonian in the form 



where Tat = ^[Fa,r{,] is used. 



n^=b±ip)-T, 



(13) 
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with 



b±ip) = 



VSpxPy 

Tvpz 



(14) 



From this, one can define the U{1) Berry flux, which is 
analogous to a magnetic field in momentum space. 



6^ 



1 

8^ 



(15) 



where b± — b±/\b±\, and the derivatives are with respect 
to p. By construction, the magnetic field is divergence- 
less, d^B'^ — 0, away from points of singularity where 
h± vanishes. However, the band touching point is such 
a singularity, and it is in fact a source of Berry flux. To 
see this, one may compute the integral of the flux of 
through a sphere around one of the touching points. It 
is straightforward to compute this integral, and by doing 
so one finds 



d^B^ip) = ±25[p). 



(16) 



Thus each band touching is a source of two quanta of 
Berry flux, and can be therefore considered a double Weyl 
point. The net Berry flux of both double Weyl points 
added together vanishes, which is required as the Bril- 
louin zone is a closed manifold without any boundary 
through which a net flux may escape. 

The Berry flux of each double Weyl point is directly 
related to the Hall conductivity. Indeed, one can show 
quite generally that the Hall conductivity is proportional 
to separation of the points, 28, 
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where the first factor of 2 is due to the doubled nature 
of the Weyl points. This is valid at zero temperature 
when the Fermi level is at the energy of the double Weyl 
points. Neglecting the small corrections to the exponent 
due to Coulomb interactions, we have K ^ \fmH , which 
explains the estimate in Eq. ([t]) of the main text. 



RG equations 

Here we describe the renormalization of the fermionic 
propagator, which was not presented in the main text. 
We start with the Hamiltonian 



Ho(fc) 



2M( 



-r" + e,(fc)r'^. 



where we defined r° as the 4x4 unit matrix to empha- 
size the matrix structure. The corresponding effective 



Lagrange density is 

CcS = -0^(-*'^rir" + Hoik) - T,f{k,iujn))ip- 

Here we have included the one-loop self energy needed 
for the £ expansion. It is given by 



^^{k,iujn) 



(17) 



' I TY,Gf{k + q, 



where G'7^(fc,iw) = -iujT^ +Ho{k), G~\q,in) = cqq^ 
are used. Note that the bosonic propagator does not 
have frequency dependence since it represents the instan- 
taneous Coulomb interaction. By integrating out high 
momentum degrees of freedom, we investigate the renor- 
malization of fermion parameters. 

Frequency and momentum dependences of the self en- 
ergy determine anomalous dimensions of the fermion 
field and mass terms, respectively. At one-loop level, 
the fermion self energy does not have frequency depen- 
dence due to the instantaneous Coulomb interaction, so 
no anomalous dimension of the fermion field appears. 
Moreover, the frequency-independence makes one mass 
coupling (1/A/o) unchanged at one loop level and it is be- 
cause frequency and this mass term have the same unity 
matrix structure. Note that both corrections naturally 
appear in higher-loop contributions. 

On the other hand, non-trivial renormalizations of 
isotropic mass (m) and anisotropic mass {Mc) appear 
even at one-loop level. For example, the renormalization 
of isotropic mass term can be read off by calculating the 
third component (F^) of the self energy. 
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It can be understood as the self energy corrrection to the 
original Hamiltonian, 



e3(fc)^e3(fc)-7'(0)fc.fcy 

\ Mc p 



(18) 



In the last line, the momentum integration is done for 
fi < q < A. The correction to the isotropic mass term, 
S(l/m) — ^/i(-^) log(^) can be regularized in the stan- 
dard way. Along the same line, one can determine the 
renormalization of the anisotropic mass term by calcu- 
lating the fourth component (F"*) of the self-energy. 

These considerations lead to the full RG equations as 



7 



follows. 



dVm' 

d . m ^ 



, m \ , 1 , 
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Here the following functions of anisotropy are introduced. 
d^k- 



/i(^) - 
h{x) - 

h{^) - 

Fi{x) 
F2{x) 



1 



1 
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3 

8^ 



dn,. 



dn,- 



- x{2 + x)pc{k) 



I + x(2 + x)pc{k) 

^ih.2(i.2 _ i.2\2 



I + x{2 + x)pc{k) 



2/3(2^) 



N, 



{l + x){h{x)~h{x)) , 



where Pc{k) = J2'i=i — ^fk^ is used. The first 

two functions, /i and /2, are from the fermion self energy 
and the third function, /a, is from the boson self-energy 
calculation. The asymptotic form of /s is shown above 
while its complete form can be obtained by evaluating the 
bosonic self energy function directly. Numerical evalua- 
tions of these functions are illustrated in Figs. [2]and|3] 
Useful numerical values are 
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/l(0)=/2(0) 

lim fi{x) = 

lim Mx) ~ 
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lim f2{x) 



F2(0) =0 
1 , 



F^iO) - 0.152 



/3(0) 




gestions as to how to apply the e expansion to Gamma 
matrices. [22] 

In this letter, we use one of possible regularization 
schemes; the angle dependence is first integrated in the 
original dimension and then dimensional regularization 
is used for the integrals over the magnitude of momen- 
tum. One reason for using this scheme is our need to keep 
track of the cubic anisotropy, which is specifically defined 

/llxj./ilxj 




FIG. 2. Numerical evaluations of the renormalization func- 
tions, fi{x), f2{x). The solid (red) and dotted (blue) lines are 
for fi{x) and /2(x) (/i(0) = ,f2(0) = 8/15). One can show 
that fi{x) — f2{x) decays as 1/x in the large x limit. 



FIG. 3. Numerical evaluations of the universal function, 
/a (a;). One can see that /a (a;) becomes saturated as x becomes 
large. 



We note that there is one subtle issue in the e expan- 
sion for regularization of the loop corrections. Since the 
gamma matrices are defined in specific dimensions, there 
is some ambiguity for the choice of gamma matrices in 
dimensional regularization. There have been some sug- 



in three spatial dimensions. Formally, the momentum 
shell integration gives X\e-di (f^ys^ Note that 

different regularization schemes only change numerical 
coefficients and do not change important physics. 



